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Abstract 

Background: The habenula and the thalamus are two critical nodes in the forebrain circuitry and they connect the 
midbrain and the cerebral cortex in vertebrates. The habenula is derived from the epithalamus and rests dorsally to 
the thalamus. Both epithalamus and thalamus arise from a single diencephalon segment called prosomere (p)2. Shh 
is expressed in the ventral midline of the neural tube and in the mid-diencephalic organizer (MDO) at the zona 
limitans intrathalamica between thalamus and prethalamus. Acting as a morphogen, Shh plays an important role in 
regulating cell proliferation and survival in the diencephalon and thalamic patterning. The molecular regulation of 
the MDO Shh expression and the potential role of Shh in development of the habenula remain largely unclear. 

Results: We show that deleting paired-box and homeobox-containing gene Pqx6 results in precocious and 
expanded expression of Shh in the prospective MDO in fish and mice, whereas gain-of-function of pax6 inhibits 
MDO shh expression in fish. Using gene expression and genetic fate mapping, we have characterized the expression 
of molecular markers that demarcate the progenitors and precursors of habenular neurons. We show that the 
thalamic domain is shifted dorsally and the epithalamus is missing in the alar plate of p2 in the Pax6 mutant mouse. 
Conversely, the epithalamus is expanded ventrally at the expense of the thalamus in mouse embryos with reduced 
Shh activity. Significantly, attenuating Shh signaling largely rescues the patterning of p2 and restores the 
epithalamus in Pax6 mouse mutants, suggesting that Shh acts downstream of Pax6 in controlling the formation of 
the habenula. Similar to that found in the mouse, we show that pax6 controls the formation of the epithalamus 
mostly via the regulation of MDO shh expression in zebrafish. 

Conclusions: Our findings demonstrate that Pax6 has an evolutionarily conserved function in establishing the 
temporospatial expression of Shh in the MDO in vertebrates. Furthermore, Shh mediates Pax6 function in regulating 
the partition of the p2 domain into the epithalamus and thalamus. 

Keywords: Organizer activity. Signaling, Forebrain, Habenula, Thalamus, Mice, Zebrafish 



Background 

The habenula, which is a paired midline structure resid- 
ing at the posterior and dorsal surface of the thalamus, 
is present in virtually all vertebrates. The habenula re- 
ceives inputs from the forebrain limbic system and the 
basal ganglia through the stria medullaris and projects 
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to the monoaminergic nuclei in the midbrain and hind- 
brain via the fasciculus retroflexus (FR) [1,2]. The neural 
connections of the habenula imply that it plays an im- 
portant role in modulating emotion, motivation and re- 
ward values. Indeed, dysfunction of the habenula has 
been implicated in psychiatric disorders, such as depres- 
sion, schizophrenia and drug-induced psychosis [2,3]. 
Despite the importance of the habenula, relatively little 
is known about the molecular control of the specifica- 
tion and differentiation of this structure in vertebrates. 

Based on gene expression patterns and morphological 
landmarks, the caudal forebrain is divided into three 
segments, prosomeres (p)l, 2 and 3 [4,5]. The pi and p3 
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domains generate the pretectum and prethalamus, re- 
spectively, while the alar plate of p2 produces the epi- 
thalamus and the thalamus. The epithalamus gives rise 
to the habenula, together with the pineal gland and 
choroid plexus. Using genetic fate mapping, we have 
recently shown that thalamic neuron precursor cells 
express homeobox gene Gbx2, and the Gbx2 lineage 
defines a sharp compartment boundary between the 
habenula and thalamus in the mouse embryo [6], dem- 
onstrating a clear segregation of habenular and thalamic 
neurons. In parallel, we have demonstrated that the 
segregation between pi and p2 cells is regulated by the 
cell adhesion factor protocadherin 10b in zebrafish [7]. 
However, how the epithalamus diverges from the thal- 
amus and pretectum is largely unknown. 

Sonic hedgehog (Shh) is expressed in the ventral mid- 
line of the neural tube and functions as a morphogen to 
control the dorsoventral patterning of the entire central 
nervous system [8,9]. Uniquely, Shh is expressed in a 
transverse band at the border between p2 and p3, called 
the mid-diencephalic organizer (MDO) located at the 
prospective zona limitans intrathalamica (ZLI) [10,11]. 
Therefore, the developing diencephalon receives Shh sig- 
nals from both dorsoventral and anteroposterior directions. 
Experiments in chicks and zebrafish have demonstrated 
that Shh is the key component of the MDO that controls 
development of the thalamus and the prethalamus in verte- 
brates [12-16]. Genetic studies in mice and fish have re- 
cently shown that Shh signaling regulates the division of 
thalamic progenitor regions into rostral (rTh) and caudal 
thalamic (cTh) domains [17-20]. Enhancing Shh activity 
shifts the border between rTh and cTh caudally, whereas 
decreasing Shh signaling shifts the border anteriorly 
[18,21,22]. Interestingly, cells in the caudal area of the 
thalamus seem refractory to alterations of Shh activity 
[18], and the effect of changed Shh function on epitha- 
lamic development has not been characterized in detail in 
the previous studies. 

Pax6, which encodes a transcription factor containing a 
paired domain and a homeodomain, is broadly expressed 
in the developing forebrain of vertebrates [23-25]. In the 
absence of Pax6y although the regionalization of the 
diencephalon is initially preserved, formation of the pro- 
somere boundary and axonal projections from the 
thalamus are severely disrupted in mice [25-28]. The 
formation of the epithalamus has never been charac- 
terized in Pax6 'dQiiciQnt vertebrates. However, it has 
been reported that the pineal gland is missing or 
greatly reduced in people lacking a functional copy of 
PAX6 [29,30], suggesting that Pax6 protein may play a 
crucial role in the development of the epithalamus in- 
cluding the habenula. Interestingly, Shh expression and 
perhaps its activity are enhanced at the MDO of Pax6' 
deficient mouse embryos at embryonic day ((E) 10.5) 



[26,28]. Antagonistic interactions between Pax6 and 
Shh are known to play a crucial role in patterning the 
spinal cord and the telencephalon [8,31]. However, 
it remains to be explored how Pax6-Shh interaction 
regulates the development of the mid-diencephalon, 
particularly the habenula. 

The goal of this study is to investigate the molecular 
control of habenula development. To this end, we first 
analyzed the expression patterns of molecules that de- 
fine different populations of habenular cells. Using these 
newly discovered habenular markers, we examined the 
genetic interactions between Pax6 and Shh in the con- 
trol of habenular development. We demonstrate that 
Pax6 regulates the temporospatial expression of Shh at 
the MDO. Furthermore, Shh signaling plays a crucial 
role in regulating the specification of the habenula. Our 
results provide novel information on the development of 
the habenula and the MDO. 



Results 

Characterization of progenitors and precursors of 
habenular neurons 

To investigate the developmental program that controls 
the subdivision of p2 into the epithalamus and the thal- 
amus, we began by looking for molecular markers that 
define the epithalamus, particularly the habenula, in 
mouse embryos. We performed a detailed analysis of a 
variety of known developmental regulators that are 
differentially expressed in the diencephalon at E12.5, a 
stage when most habenular and thalamic neurons are 
being generated. As described previously [17], proneuro- 
nal gene Neurog2 (previously named Ngn2) is broadly 
expressed in progenitors located in the thalamic ven- 
tricular zone (VZ) and also basal progenitors in the cTh 
(Figure lA). Neurog2 transcripts are also present in the VZ 
dorsal to the cTh, suggesting that Neurog2 is also expressed 
in epithalamic progenitors (Figure lA). Irxl, which encodes 
a protein belonging to the Iroquois homeodomain family, is 
expressed in the dorsal side of the diencephalon [32]. At 
E12.5, transcripts of Irxl, as well as those of paralogous 
genes Irx2 and IrxS, are abundantly present throughout the 
epithalamus, pretectum and prethalamus, while a lower 
level of expression is seen in the VZ of the cTh (Figure IB 
and data not shown). Homeobox gene Dbxl is initially 
expressed in p2 in a dorsal^^^^-ventral^^^ gradient [18,33]. 
By E12.5, Dbxl transcripts are mainly restricted to the VZ 
of the presumptive epithalamus and are barely detectable 
in the cTh (Figure IC). In addition to its strong expression 
in the prethalamus, a weaker expression domain of Wnt7b 
is present in the VZ of the epithalamus overlapping with 
that of Dbxl (Figure ID). Collectively, our data suggest that 
the habenular progenitor domain is marked by the com- 
binatorial expression of Neurog2, Irxl, Dbxl and Wnt7b, 
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Figure 1 Identification of molecular markers for the habenular precursor domain. (A-l) In situ hybridization (ISH) for tine indicated marl<ers on 
coronal sections of wildtype embryos at embryonic day (E)12.5. Tlie inset in I sliows ISH for Pax3. Brackets indicate the presumptive habenular domain; 
the dashed line in I demarcates the border between the pretectum and epithalamus. (J,K) Schematic representation of coronal and sagittal sections of 
El 2.5 brains, respectively. The dashed line in K indicates the plan of sections. (L-M) Contribution of descendants (RFP^) of ^t\// -expressing cells that are 
labeled at El 2.5 to the medial habenula (MHb) in Etvl^^^''^^^; R26R^^^^^ at post-natal day (P)0 (L) or at P22 (M). The distribution of the fate-mapped cells is 
summarized in L' and the MHb is marked by Pou4fl immunoreactivity (M). cTh, caudal thalamus; ^h, epithalamus; LHb, lateral habenula; PrT, pretectum; 
rTh, rostral thalamus; ZU, zona limitans intrathalamica. Scale bar in (A): 200 \im (A-l); 250 |jm (L); 196 pm (M). 
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We next analyzed genes that mark postmitotic habe- 
nular neurons. As described previously, Pou4fl (also 
known as BrnSa) is specifically expressed in the habe- 
nula [33] (Figure IE). Axonal guidance genes RoboS and 
Neuropilin 2 {Nrp2) are expressed in the habenula 
[33,34]. Transcripts of RoboS and Nrp2 were detected in 
cells that emerged from the habenular progenitor do- 
main marked by Dbxl and Wnt7b (Figure IC-G). Etvl^ 
which encodes a member of the ETS family of transcrip- 
tion factors, is expressed in the prospective habenula at 
E14.5 [33]. At E12.5, Etvl transcripts were detected in 
cells residing between the lateral-most area that was posi- 
tive for Pou4fh RoboS and Nrp2 and the Dbxl^lWntVb^ 
VZ (Figure IH). Although it has been suggested that 
Bhlhe2S (previously called Bhlhb4) is expressed in the area 
caudal and dorsal to the thalamus including the epithal- 
amus [35], Bhlhe2S transcripts were mainly found in the 
lateral wall of the pretectum residing in an area dorsal and 
caudal to the Pou4fl-\- domain (Figure 1E,I). Furthermore, 
the Bhlhe2S^ domain was positioned parallel but more lat- 
eral to the expression domain of PaxS, which demarcates 
the pretectal VZ (Figure II and inset) [36,37]. 

To forge the link between the gene expression pattern 
in the embryo and the mature habenular nuclei after birth, 
we performed genetic fate mapping to examine the devel- 
opmental fate of £^vi -expressing cells in the mouse em- 
bryo. We crossed the Etvf^^^^ knock-in mouse strain 
(JAX013048) and a cre-reporter line, R26lf^^, in which 
robust red fluorescence protein (RFP) will be permanently 
expressed upon cre-mediated excision of a STOP cassette 
flanked by two loxP sequences upstream of a cDNA 
encoding tdTomato RFP [38]. No RFP-labeled cells were 
detected in Etvr'^"^^^; R26]f'^''^^ embryos that did not 
receive tamoxifen (data not shown), demonstrating that 
creER-mediated recombination is dependent on tamoxifen 
administration. To label ^^vi -expressing cells, we ad- 
ministered tamoxifen to time-pregnant females carrying 
^f-yjcreERU, j^26R^^^'^^ embryos, and examined the dis- 
tribution of RFP-labeled cells at postnatal day (P) 0 and 
22. In the caudal forebrain, descendants of £^vi -ex- 
pressing cells labeled at E12.5 were mostly found in 
the medial habenular nuclei, which were delineated by 
Pou4fl immunoreactivity, in Etvl'''^^^^; R26R^^^ mice 
at PO and P22 (Figure IL-M). Therefore, the Etvl ex- 
pression domain in the diencephalon at E12.5 identifies 
precursors of the medial habenular nuclei. 

In summary, we have identified several molecular 
markers that are expressed in distinct but overlapping 
domains of the prospective habenula at E12.5 (Figure IJ). 

Loss of Pax6 results in expansion of the pretectum and 
thalamus at the expense of the epithalamus 

To investigate if Pax6 is essential for habenular develop- 
ment, we characterized the diencephalic phenotype in 



mouse embryos homozygous for a spontaneous Small 
eye {Sey) mutation, which is a null mutation in the Pax6 
gene [39]. The transcripts of Dbxl and Wnt7b were 
completely missing in p2 oi Paxe^'^^^'y embryos at E12.5 
(Figure 2H,I). Furthermore, Neurog2 expression in the 
presumptive habenular progenitor domain was also ab- 
sent (Figure 2J). In concurrence with the abnormal specifi- 
cation of the habenular progenitor domain, transcripts of 
Irxl, Pou4fl, Nrp2 and Etvl were mostly absent in the pro- 
spective habenular precursors of Pax6^^^^^^^ embryos at 
E14.5 (Figure 3B,F,J,N). These observations indicate that 
loss of Pax6 disrupts the specification and differentiation of 
the habenula. 

We next examined how loss of Pax6 altered the pat- 
terning of the diencephalon, particularly the formation 
of the thalamus and the pretectum, as they reside in the 
ventral and caudal sides of the habenula, respectively. 
OligS is normally expressed in the thalamus (rTh and 
cTh) and the ZLI [17] (Figure 2D). In Pax6^'^^^'^ em- 
bryos at E12.5, OligS expression was maintained and its ex- 
pression domain appeared expanded dorsally (Figure 2K). 
Because of the flexure of the neural tube, the pretectum, 
which is demarcated by the expression of PaxS and Ascll, 
appeared dorsal to the epithalamus on a coronal section 
(Figure 2E,F). Immunofluorescence for Dbxl and Pax3 
showed that these two molecules are expressed in two jux- 
taposed domains, corresponding to the habenula and the 
pretectum, respectively, at E12.5 (Figure 2G). In Paxe^^^''^^^ 
embryos, the expression domains of PaxS and Ascll were 
expanded anteriorly/ventrally so that the Ascll and PaxS 
expression domain abnormally opposed or partially 
overlapped with that of the OligS^ domain at E12.5 
(Figure 2K-M). Remarkably, probably because of the 
protein perdurance, Dbxl protein was still detectable 
in the presumptive habenula of Pax6^^^^^^^ embryos 
despite the lack of Dbxl transcripts by E12.5 (Figure 2H,N). 
However, many of these Dbxl^ cells abnormally expressed 
Pax3 (Figure 2N). This suggests that some habenular pro- 
genitors may be initially formed but adopt a pretectal fate 
in the absence of Pax6, Collectively, our data suggest that 
the thalamic and pretectal progenitor domains are respect- 
ively expanded dorsally and anteriorly at the expense of the 
epithalamus in the absence of Pax6, 

To confirm the abnormal expansion of thalamic and 
pretectal domains in Pax6^^^^^^^ embryos, we analyzed 
markers for thalamic and pretectal neurons. Homeobox 
gene Gbx2 is broadly expressed in postmitotic thalamic 
neurons and the G^:v2-lineage exclusively contributes to 
cTh-derived nuclei [6,40]. Taking advantage of enhanced 
green fluorescence protein (EGFP) expression from the 
Gbx2^''^^^ locus, which contains a creER-ires-EGFP cas- 
sette in the 5 '-untranslated region (UTR) of the Gbx2 
gene [6], we examined the distribution of thalamic neu- 
rons m Pax6^'y^^'y embryos at E14.5. Double labeling of 
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Figure 2 Pattern formation of the diencephalon in embryos without Pax6 or with reduced Shh signaling. (A-F, H-M, and 0-T) ISH on 

adjacent coronal sections of El 2.5 embryos. Probes and genotypes are indicated on tine top and left of the panel. (G, N, and U) 
Immunofluorescence for Dbxl and PaxS on coronal sections of El 2.5 embryos. Red dashed lines in C and D demarcate the border between the 
epithalamus and the thalamus; dashed line arrows indicate the expansion of the particular progenitor domain; black arrowheads and brackets in 
C, J, Q and M show the normal and the expanded ZLI, respectively; black arrows in F and T indicate the rTh domain; red brackets demarcate the 
prospective habenular domain; white arrowheads in N show cells positive for both Dbxl and PaxS. Note that the apparent lack of transcripts of 
Pqx3 and Ascl in the pretectum in S and T is due to an anteriorly tilted sectioning plane. Scale bar in (A): 200 pm (A-F, H-M, and 0-T); 241 pm 
(G, N, and U). E, embryonic day; ISH, in situ hybridization; rTh, rostral thalamic; ZLI, zona limitans intrathalamica. 



GFP and neurofilament (NF) showed that the posterior 
border of the GFP^ thalamus was marked by the NF^ FR 
tract, which is formed by the habenular efferents [41] 
(Figure 4A,A')- Pretectal neurons are mostly gamma- 
aminobutyric acid (GABA)ergic and extend axonal fibers 
to the dorsal midline forming the posterior commissure 
(Figure 4A). In Paxe^'^^^'^ embryos at E14.5, GFP"" cells 
encroached the presumptive habenula and were inter- 
mingled with GABA^ neurons (Figure 4B,B')- The FR tract 
was replaced by defasciculating axonal fibers running from 
the GFP"^ domain to the ventral midbrain (Figure 4B,B')« 
Interestingly, these axonal fibers were mostly positive for 
GABA, resembling posterior commissural axons rather 
than the FR tract (Figure 4B). Collectively, these data 
suggest that thalamic and pretectal neurons occupy the 
presumptive habenula in Pax6^^^^^^^ mutants. 



In summary, our findings demonstrate that Pax6 is 
essential for the development of the habenula. In the 
absence of Pax6y the habenula is lost in association with 
the dorsal expansion of the thalamus and anterior ex- 
pansion of the pretectum. 

Pax6 is required for the establishment of the aniage of 
the MDO 

We next investigated if Pax6 has a conserved function 
within the vertebrate lineage. Therefore, we analyzed 
Pax6 function during mid-diencephalic development in 
zebrafish. First, we analyzed the expression pattern of a pre- 
thalamic marker, dlx2a, and pretectal marker, proxl [16]. 
We found that both expression domains expand towards 
the MDO/thalamus aniage and simultaneously the MDO/ 
thalamus primordia shrink in embryos overexpressed with 
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Control P^^^^^^f^l Shh^'^ <^ PsxB^^^y; Shh^ 




Figure 3 Opposing function of Shh and Pax6 in regulating formation of tlie liabenula. (A-P) ISH for Irxl, Pou4fl, Nrp2 and Etvl on sagittal 
(A-L) and coronal (M-P) sections of El 4.5 embryos of indicated genotypes. Insets show coronal sections of El 4.5 embryos. Red dashed lines 
indicate the posterior border of the thalamus demarcated by the PR tract (arrow); green dashed lines demarcate the posterior limit of the 
habenula; red brackets show the habenular nuclei; red and green asterisks indicate the loss and ectopic expression of habenular markers, 
respectively; empty arrowheads show abnormal defasciculating axonal fibers replacing the PR tract. Scale bar in (A): 200 pm (A-L); 170 pm (M-P); 
71 1 |jm (insets). E, embryonic day; PR, fasciculus retroflexus; ISH, in situ hybridization. 



Pax6 mRNA (n = 7/11; Figure 5A,C,G). To block Pax6 
function in zebrafish, we performed a morpholino 
oligomer-based double knock-down approach for both 
Pax6 homologs, Pax6a and Pax6b, which are expressed 
early in the diencephalon of zebrafish embryos [42]. 
Consistently, knock-down of Pax6 leads to a broaden- 
ing of the MDO/thalamus anlage (n = 6/10) and in a 
few cases the pretectal expression of proxl becomes 
down-regulated (n = 3/10; Figure 5E), Next, we ad- 
dressed the expression of otx2, a marker of the primordium 
of the MDO/thalamus, and irxlk which is expressed in the 
thalamus (Figure 5B,G). We found that irxlb shifts anteri- 
orly into the territory of the MDO anlage in embryos in 
which pax6 mRNA is overexpressed (n = 8/13; Figure 5D), 



suggesting that Pax6 represses the MDO primordium. 
Consistently, knock-down of Pax6 function leads to a 
broadening of the MDO primordium, whereas the thal- 
amus seems to be less affected (n = 6/13; Figure 5F). This 
suggests that Pax6 restricts the anteroposterior width of the 
MDO primordium in zebrafish. 

Pax6 is essential for regulating the expression of Shh in 
the prospective MDO 

Pax6 is broadly expressed in the diencephalon before 
E10.5 and maintained in the habenula and prethalamus 
at E12.5 [23]. Previous studies have shown that loss of 
Pax6 leads to expanded expression of Shh in the MDO 
[26,28], suggesting that Pax6 plays an important role 
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Figure 4 The partition of p2 to form thie hiabenula and thie thalamus is controlled by the interaction between Pax6 and Shh. (A-C) 

Immunofluorescence for GABA, GFP, and neurofilament (NF) on sagittal sections of El 4.5 embryos of indicated genotypes. (A'-C) the blue 
channel of images A-C to show for NF staining. The dashed line outlines the habenula (Hb); arrows indicate the FR tract; the asterisk shows 
defasciculating axonal fibers in place of the FR tract. Note that the FR tract becomes enlarged and supernumerary in Shh^^^^^ embryos 
(see inset in C). Scale bar in (A): 200 pm (A-C). E, embryonic day; FR, fasciculus retroflexus; GABA, gamma-aminobutyric acid; GFP, green 
fluorescent protein. 



in regulating the MDO Shh expression. We confirmed 
that the Shh expression domain in the MDO was 
indeed expanded in Pax6^^^^^^^ embryos at El 2. 5 
(Figure 6A,B). Furthermore, the expression domain of 
Ptchl, which is considered as a readout gene of Shh 
activity [43], was also expanded in both the thalamus 
and prethalamus in Pax6^^^^^^^ embryos at El 2. 5 
(Figure 6E,F), demonstrating enhanced Shh signaling 
in the diencephalon lacking Pax6, It has been shown 
that activating Shh by forced expression of a constitu- 
tively active Shh receptor Smoothened (Smo) through- 
out the neural tube causes expansion of the rTh at the 
expense of the cTh [18]. In accordance with these findings, 
we found that the ZLI and the rTh, which are demarcated 
by the adjacent AsclF /Neurog2^ and Ascll^/Neurog2~ do- 
mains, respectively [18], were noticeably enlarged in 
Pax6''y^''y embryos at E12.5 (Figure 2C,F,J,M). Altogether, 
these results demonstrate that the loss of Pax6 results in 
an expanded Shh expression domain and enhanced Shh 
activity at the ZLI. 

We next investigated whether Pax6 was required for 
the establishment of Shh expression in the ZLI. Shh is 
broadly expressed in the basal plate of the future 
diencephalon as early as E8.5 [20]. A wedge-shaped 
Shh^ domain corresponding to the prospective ZLI was 
detected in wildtype embryos at E10.5 but not at E9.5 
(Figure 61 and data not shown). By contrast, the MDO 
Shh expression domain was clearly visible in Pax6^^^^ 
embryos by E9.5 (n = 3/3, Figure 6J), demonstrating 
that the loss of Pax6 results in precocious induction of 
MDO Shh expression in mice. 



Loss of Shh activity leads to expansion of the epithalamus at 
the expense of the thalamus at El 2.5 

To investigate whether the loss of the habenula could be 
attributed to altered Shh function in Pax6^^^^^^^ embryos, 
we first sought to determine if Shh is essential for habe- 
nular development. There is decreased cell proliferation 
and increased cell death in the entire diencephalon in 
mouse embryos homozygous for a Shh-miW mutation at 
E9.5 [44] . To investigate the patterning role of Shh with- 
out the complication of growth defects in the dienceph- 
alon lacking Shh before E10.5, we performed analyses in 
mice homozygous for a Shh knock-in allele (designated 
as Shh^^), which produces Shh::GFP fusion protein with 
reduced Shh activity [45]. In Shh^^^^^ embryos at E12.5, 
although Shh expression was unchanged in the basal 
plate and ZLI, the expression domain of Ptchl was no- 
ticeably reduced in the thalamus at E12.5 (Figure 6C,G), 
indicating attenuated Shh activity in the p2 area in 
Shh^^^^^ embryos. Remarkably, the expression domain 
of Dbxl and Wnt7b in the presumptive epithalamus was 
expanded ventrally and overlapped with that of Neurog2 
and OligS at E12.5 (Figure 20-R). At E14.5, the expres- 
sion domain of Irxl, Pou4fl, Nrp2 and Etvl was notice- 
ably enlarged at the expense of the thalamus in Shh^^^^^ 
embryos (Figure 3C,G,K,0). Immunofluorescence for 
GFP and NF in Gbx2'''^^^^; Shh^^^^^ embryos at E15.5 
revealed that the GFP^ thalamus was evidently reduced 
in size compared with that of the control (Figure 4C). 
Furthermore, the habenula was expanded ventrally, and 
noticeably enlarged and multiple FR tracts were often 
detected in Shh^^^^^ embryos (Figure 4C,C' and the 
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Figure 5 Pax6 is essential for thie establishiment of thie mid-diencephialon in zebrafisii. (A-F) Whole mount double in situ hybridization 
(ISH) for indicated markers of 28 hours post-fertilization (hpf) zebrafish embryos of indicated experimental procedure. Double red arrows mark 
MDO thalamic primordial. In B, D and F, black dashed line marks the border between MDO and the prethalamus, and white dashed line marks 
the border between MDO and thalamus. (G) Schematic representation of a sagittal section of zebrafish forebrain at 28 hpf. cTh, caudal thalamus; 
ETh, epithalamus; MDO, mid-diencephalic organizer; PrT, pretectum; PTh, prethalamus; rTh, rostral thalamus; ZLI, zona limitans intrathalamica. 



inset). Our results demonstrate that Shh plays an important 
role in positioning the border between the epithalamus and 
the thalamus. Reducing Shh activity leads to enlargement 
of the habenula at the expense of the thalamus. 



Reducing Shh signaling partially restores the formation of 
the habenula in Pax6^^^^^^^ embryos 

To define the epistatic relationship between Shh and Pax6 
in regulating habenular development, we investigated 
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Control 




Figure 6 Loss of Pax6 results in precocious formation and expanded expression of Shh in the MDO. (A-H) ISH for Shh and Ptchl on 
coronal sections of El 2.5 embryos of the indicated genotypes. (I and J) Wliolemount ISH for Shh in control (I) and Pox6^^^^^^^ (J) embryos at E9.5. 
Arrows denote the MDO; brackets show the expanded domain of Shh in the MDO; dashed line arrow indicates the dorsally (rostrally) and 
ventrally (caudally) expanded Ptchl expression domains in F; the red arrowhead in J indicates the precocious induction of the ZLI Shh expression 
in Pax6^^^^^^^ embryos. Scale bar in (A): 200 pm (A-H). E, embryonic day; ISH, in situ hybridization; MDO, mid-diencephalic organizer; ZLI, zona 
limitans intrathalamica. 



whether reducing Shh activity could rescue the habe- 
nula in embryos lacking Pax6 by generating Pax6^^^^^^^) 
Shh^^^^^ double mutant embryos. Similar to that found in 
PaxS^'y'^^'y embryos, the MDO Shh expression domain 
was expanded in Pax6^^^^^^^; Shh^^^^^ double mutant em- 
bryos at E12,5 (Figure 6D), However, the range and the 
level of Ptchl expression in the thalamus in Paxe^'^^^^; 
Shh^^''^^ embryos were more similar to those found in 
wildtype than those in either Paxe^^^''^^^ or Shh^^^^^ single 
mutant embryos at E12.5 (Figure 6E-H). These observa- 
tions suggest that attenuated Shh protein activity may par- 
tially offset an increased Shh transcription at the MDO 
leading to relatively normal Shh signaling in the dienceph- 
alon in PaxS^'^^^'^; Shh^^^^^ embryos. Significantly, the ex- 
pression domain of Irxl, Pouf41 and Nrp2 was largely 
restored in the presumptive habenular region in Pax6^^^^^^^\ 
Shh^^^^^ double mutant embryos at E14.5 (Figure 3D,H,L, 



n >3 for each probe). These results indicate that Shh acts 
downstream of Pax6 in controlling the formation of 
the habenula. Despite the rescue of some habenular 
markers, Etvl expression and the FR tract were still absent 
in Paxe^'^^^'y; Shh^^^^^ embryos at E14.5 (Figure 3D,H,L,P). 
Therefore, in addition to its important role in restricting 
Shhy Pax6 may have an additional role in directly regulating 
habenular formation. 

Pax6 limits Shh expression in the MDO independent of 
Irx1b function in zebrafish 

Next, we investigated the role of Pax6 in the regulation 
of the MDO Shh expression in zebrafish. Overexpression 
of Pax6 led to a strong decrease of shh expression in the 
MDO in zebrafish (n = 16/26; Figure 7A-B'), consistent with 
the loss of Otx2 in the MDO primordium (Figure 5D). 
Similar to the loss of Pax6 in mice, blockage of the function 
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Figure 7 Pax6 regulates the expression domain of shh expression in the MDO in zebrafish. (A-H) Lateral views of wholemount ISH for shh 
in 28 lipf zebrafisli embryos of tine indicated experimental procedure. (A-D) red arrows mark the MDO, whereas red arrowheads mark the 
position of the reduced MDO. (A^B',C,D') show a coronal section of the prosencephalon with anterior to the bottom, hpf, hours post- 
fertilization; ISH, in situ hybridization; MDO, mid-diencephalic organizer. 



of both zebrafish homologues pax6a and pax6b by mor- 
pholino oligomers led to an expansion of the MDO shh 
expression domain (n = 6/12; Figure 7C-C'); however, a pre- 
cocious formation of the shh positive organizer was not ob- 
served in fish (data not shown). Co-injection of human 
PAX6 mRNA restored the MDO shh expression in pax6 
double morphant zebrafish embryos (n = 3/7; Figure 7D,D'), 
demonstrating the specificity of the knock-down. In fish 
and mice, it has been shown that irxl controls the posterior 
limit of the MDO shh expression domain [15,46]. We 
found that overexpression of PAX6 led to an anterior shift 
of irxlb expression into the primordium of the MDO 
(Figure 5D). Therefore, we asked whether the down- 
regulation of shh expression in Pax6-overexpressing em- 
bryos is mediated by irxlb or directly controlled by Pax6. 
Depletion of the irxlb function in zebrafish results in ex- 
panded shh expression in the MDO, similar to that found 
in pax6 double morphant embryos (n = 6/9; Figure 7C,G). 
To determine the epistatic relationship between pax6 and 



irxlb in regulating shh expression, we injected pax6 
mRNA in irxlb-MO embryos and found that the MDO 
shh expression domain was greatly reduced, similar to that 
found in paxS-overexpressing embryos (n = 4/8; Figure 7H). 
These data show that pax6 acts downstream or in parallel 
with irxlb in restricting the MDO shh expression in fish. 

Collectively, our results demonstrate that Pax6 plays 
an evolutionarily conserved function in establishing the 
temporal and spatial expression domain of Shh in the 
prospective MDO organizer. 

Shh function from the MDO influences habenula 
formation 

In the next set of experiments, we investigated if Pax6 
and Shh regulate the formation of the epithalamus in 
zebrafish in a way similar to that found in the mouse 
embryo. We used otxS as an early marker for the pineal 
gland, a key derivative of the epithalamus [47] and 
cxcr4b as a marker for habenula precursors [48]. We 



Chatterjee et at. BMC Biology 2014, 12:13 
httpy/www.biomedcentral.com/l 741 -7007/1 2/1 3 



Page 11 of 15 



found that double knock-down of pax6a/pax6b led to a 
reduced expression of otxS (n = 7/15; Figure 8A-B') and 
cxcr4b (n = 11/19; Figure 8G,H'). Conversely, overexpres- 
sion of pax6 mRNA led to an increase in the expression 
domain of otxS (n = 4/7; Figure 8C,C') and cxcr4b (n = 61 
10; Figure 81,1')- This indicates that Pax6 positively regu- 
lates the size of the epithalamus in fish similar to that 
found in mouse embryos. We reasoned that the influ- 
ence of Pax6 on epithalamus development could be me- 
diated by Shh signaling from the MDO as found in 
the mouse embryo. Therefore, we analyzed epithalamus 
formation in embryos carrying a mutation in the Shh 
co-receptor smoothened (referred to as Smu^^^^^^^^ mu- 
tants) [49]. We found an increase in the epithalamic ex- 
pression of OtxS in Smu^^^^^^^^ mutants (Figure 8D,DM,J'). 
Next, we blocked Pax6 function in the Smu^^^^^^^^ mu- 
tant background and found that in the absence of Shh sig- 
naling Pax6 knock-down does not alter the size of the 
epithalamic otxS expression (Figure 8E,E') and cxcr4b ex- 
pression domain (Figure 8K,IC) compared to Smu^^^^^^^ 
mutant embryos (Figure 8D,D',J,J'). Similarly, we found that 
overexpression of Pax6 in Smu^^^^^^ did not change the 
epithalamus phenotype (Figure 8F,F',L,L') compared to 
Smuf"^^^^^^ mutant fish. This suggests that Pax6 influences 
the development of the epithalamus - here the pineal 



complex and habenula - via regulation of Shh signaling at 
the ZLI organizer in the mid-diencephalon. 

Discussion 

There is a surge of interest in the habenula because of 
its important role in cognitive and emotive behaviors 
[2,3,50,51]. However, the molecular underpinnings of the 
development of the habenula in vertebrates are still 
poorly understood. Expression profiling of postmitotic 
habenular neurons has recently demonstrated that the 
habenula represents a unique molecular territory in the 
central nervous system and is composed of heteroge- 
neous ceU types [33]. As a first step in understanding 
the specification and differentiation of habenular neu- 
rons, we have characterized the molecular markers that 
define habenular progenitors and precursors. Using gen- 
etic fate mapping, we have established a definite link be- 
tween the medial habenular nuclei and £^vi -expressing 
progenitors in the epithalamus. 

Transplantation and genetic manipulation experiments 
have demonstrated that Shh is the key molecule for 
MDO activity [12-14,16]. Experiments in the zebrafish 
have suggested that transcription factors, such otxl/otx2, 
fezf2 and irxlb, form an interactive network to define 
the competent domain for the induction of shh in the 
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Figure 8 Pax6 controls development of the epithalamus via regulation of MDO Shh function in zebrafish. (A-L) Lateral and dorsal views 
of wholemount ISH for otx5 (A-F') and cxcr4b (G-L') in 28 hpf zebrafish embryos of the indicated experimental procedure. Red arrows marks the 
cxcr4b positive habenula precursor cells, hpf, hours post-fertilization; ISH, in situ hybridization; MDO, mid-diencephalic organizer. 
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prospective MDO [16]. This model, however, is not 
completely compatible with expression data and loss-of- 
function studies in mice [11]. The molecular mechanism 
underlying the formation of the MDO Shh expression in 
mammals remains to be determined. We found that loss 
of Pax6 causes not only expanded Shh expression, but 
also precocious formation of the wedge-shaped Shh ex- 
pressing domain in the prospective ZLI in mouse em- 
bryos by E9.5, at least 24 hours before the expression 
normally commences. To the best of our knowledge, this 
is the first report of accelerated formation of the MDO 
Shh expression domain caused by a mutation in verte- 
brates. It has been shown that Shh function in the ven- 
tral region is not essential for the induction of the MDO 
Shh expression [16,20]. In line with these conclusions, 
we found that the MDO Shh expression domain was 
similarly expanded in Paxe^'^^^'^ and Paxe^'^^^'^) Shh^"^^^"^ 
double mutants. These findings demonstrate that the tem- 
porospatial expression of Shh in the prospective ZLI is in- 
dependent of Shh signals from the basal plate, but 
directly regulated by Pax6. Gain-of-function of Pax6 
results in reduced expression of the MDO shh expres- 
sion in fish (Figures 5 and 7) and chick embryos as 
described recently [22]. Therefore, Pax6 plays an evo- 
lutionarily conserved function in the establishment of 
the MDO Shh expression. 

In thalamic explants, different concentrations of Shh 
proteins induced differential expression of rTh- and 
cTh-specific markers, suggesting that Shh acts as a mor- 
phogen to pattern p2 [12]. Indeed, increasing or redu- 
cing Shh activity alters the position of the border 
between rTh and cTh in relation to the ZLI [18,20]. Al- 
though the formation of the habenula was not explicitly 
characterized, altered expression of Dbxl and OligS, in- 
dicative of abnormal p2 subdivision, was reported in 
mouse embryos with excess or reduction of Shh activity 
[18]. By analyzing mouse embryos homozygous for Sey 
or a 5'/z/z-hypomorphic mutation and Smu^^^^^^^^^ mu- 
tant embryos or fish embryos knocked-down for Pax6 
function, we showed that Shh plays a crucial role in po- 
sitioning the border between the epithalamus and cTh, 
demonstrating that Shh has a long-range effect in the 
development of the entire p2. The far-reaching activity 
of Shh may be achieved by the combined action of Shh 
signals emanated from both the basal plate and the ZLI 
cells (Figure 9A). Future studies are required to deter- 
mine if Shh signaling positions the border between the 
cTh and the habenula by directly inhibiting habenula 
fate or indirectly via promoting the cTh identity. 

Antagonizing interactions between Pax6 and Shh play 
an important role in dorsoventral patterning of the 
neural tube. In the spinal cord, loss of Pax6 leads to ex- 
pansion of the Shh expression domain in floor plate [31]. 
Furthermore, Pax6 regulates the competence of the 



tissue in response to Shh signals [22,31]. Similarly, re- 
moval of Pax6 partially rescues the medial ganglionic 
eminence marker in the forebrain of Shh mutants [8] . In 
the current study, we demonstrate a similar antagonistic 
interaction between Pax6 and Shh in patterning the p2 
alar plate. In the absence of Pax6, the development of 
habenula was blocked, while markers of the thalamus 
and pretectum were ectopically expressed in the pre- 
sumptive habenula (Figures 2, 3 and 4 and summary in 
Figure 9B). Similarly, it has been shown that the expres- 
sion domain of thalamic markers, Vmat2 and Ptn (previ- 
ously known as Hbnf), is shifted dorsally to the 
epithalamus in Pax6^^^^^^^ embryos [28]. Conversely, the 
expression of habenular markers was expanded into the 
thalamus in mouse and fish embryos with reduced Shh 
function or gain of function of Pax6 (Figures 2, 3, 8 and 
9C). In addition, we observed partial co-expression of 
pretectal and caudal thalamic markers if Pax6 is 
inactivated, which suggests that Pax6 may also be in- 
volved in demarcating the cTh/pretectum boundary. 
Importantly, attenuation of Shh partially rescued ex- 
pression of habenular markers in the presumptive epi- 
thalamus in Pax6^^^^^^^ embryos (Figures 3 and 9D). 
Furthermore, genetic experiments in zebrafish clearly 
demonstrated that shh acts downstream of pax6 in 
controlling the formation of the pineal complex and 
habenular precursors (Figure 8). Collectively, our find- 
ings demonstrate that patterning defects in the p2 do- 
main that is devoid of Pax6 can, in part, be attributed 
to the precocious formation of MDO Shh expression 
and/or enhanced Shh activity in both mouse and fish 
embryos. It is worth noting that despite the rescue of 
some habenular markers, Etvl expression was missing 
and the FR tract was also absent in Paxe^'^^^'^) Shh^^^^^ 
embryos at E14.5 (Figure 4J-L and data not shown). 
It has been shown that Etvl is a direct transcriptional 
target of Pax6 in cortical development [52]. Further- 
more, we found that deletion of Etvl resulted in abnor- 
mal formation of the medial habenular nuclei in mice 
(Guo and Li, unpublished observations). Therefore, in 
addition to its regulation of Shh, Pax6 may also have 
a direct role in regulating the differentiation of the 
habenula. 



Conclusions 

In this study, we have combined gene expression and gen- 
etic fate-mapping to define the progenitors of habenular 
nuclei. We demonstrate that Pax6 is essential for establish- 
ing the temporospatial expression of Shh in the prospective 
ZLI in vertebrates. Furthermore, we demonstrate that Shh 
has a far-reaching effect in mediating Pax6 function in con- 
trolling the p2 domain to diverge into the epithalamus, 
cTh, and rTh. 
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Figure 9 Summary of subdivision of the diencephalon through interactions between Shh and Pax6. (A-D) Schematic representations 
illustrate the partition of the diencephalon of indicated genotypes. In wildtype embryos (A), Pox6 restricts the expression domain of S/i/i in the 
MDO. Strong Shh signaling (arrows) induces rostral thalamus (rTh), while Shh appears to have a long-range effect (dashed line arrows) in inducing 
caudal thalamus (cTh) and inhibiting epithalamus (ETh). In Pqx6 loss-of-function mutants (B), the MDO Shh expression domain is enlarged in 
mouse and zebrafish embryos, and the formation of the MDO Shh expression is accelerated in the mouse. The enhanced Shh activity causes 
expanded rTh and cTh and prethalamus (PrT) at the expense of the epithalamus. In embryos with reduced Shh function or gain of function of 
Pax6 (C), the epithalamic domain is enlarged at the expense of cTh. In Pqx6 loss-of-function mutants with reduced Shh activity (D), although the 
MDO Shh expression is enlarged, the epithalamus is partially restored. PTh, prethalamus; pi -3, prosomeres 1-3. 

V ) 



Methods 

Maintenance of mouse and fish 

All animal work has been conducted according to rele- 
vant national and international guidelines. Experimental 
procedures with mouse and fish were approved by the 
Animal Care Committee at the University of Connecti- 
cut Health Center and the Regierungsprasidium Karlsruhe 
(Aktenzeichen 35-9185.64) and the Karlsruhe Institute of 
Technology (KIT). All mouse strains were maintained on a 
mixed genetic background. Noon of the day on which a 
copulatory plug was detected was designated as E0.5. Mice 
carrying the small eye mutation Pax6^^^ were identified by 
the characteristic nasal and eye defects and confirmed with 
polymerase chain reaction (PCR) analysis as described 
previously [39]. Embryos carrying the Gbx2^^^^^ and 
Shh^^^"^^^ alleles were identified by EGFP florescence in 
the spinal cord and by PCR genotyping as described 
previously [6,45]. At least three embryos of each genotype 
were examined and only the reproducible phenotype was 
described unless otherwise indicated. 

Breeding zebrafish {Danio rerio) were maintained at 
28°C on a 14 hour light /lO hour dark cycle. To prevent 
pigment formation, embryos were raised in 0.2 mM 
l-phenyl-2-thiourea (PTU, Sigma, St. Louis, MO) after 
24 hpf. The data we present in this study were acquired 
from analysis of KIT wild-type zebrafish AB2O3 as well as 



the slow muscle omitted mutant line (referred to as 
smu) carrying a mutation in smoothened homolog [49]. 

Functional analysis 

Transient knock-down of gene expression in zebrafish 
was performed as described previously [7]. We used the 
following morpholino-antisense oligomers (Gene Tools, 
Philomath, OR) at a concentration of 0.3 mM each: pax6a 
MO (5'-TTTGTATCCTCGCTGAAGTTCTTCG-3') and 
pax6b MO (5'-CTGAGCCCTTCCGAGCAAAACAGTG- 
3') [53]. MO oligomers were injected into the yolk cell 
close to blastomeres at the one-cell or two-cell stage. 

For mis-expression experiments full-length pax6 and 
irxlb was cloned into a pCS2+ vector [15] and from 
this template mRNA was synthesized in vitro (Message 
Machine Kit, Ambion, Amersham, UK.). Together with 
rhodamine dextran (MiniRuby, Invitrogen, Carlsbad, 
CA) as lineage tracer, 350 pg mRNA per embryo was 
injected into the one-cell stage. 

Immunohistochemistry and in situ hybridization 

Fish embryos were fixed in 4% paraformaldehyde/phos- 
phate-buffered saline (PBS) at 4°C overnight for further 
analysis. Whole-mount mRNA ISHs were performed as 
described previously [54]. Embryonic mouse brains were 
dissected in cold PBS and fixed in 4% paraformaldehyde 
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overnight for immunostaining or RNA ISH. Brains 
were cryoprotected and embedded in optimal cutting 
temperature (OCT) compound (TissueTek, Torrance, 
CA, USA) and sectioned in 20 (im thickness. Immuno- 
histochemistry and ISH were performed as described 
previously [55]. Detailed protocols are available on the 
Li lab website [56]. The following antibodies were used 
in this study: rabbit anti-GFP (Invitrogen, Carlsbad, CA, 
USA); mouse anti-neurofilament, anti-Pax3, and anti- 
Pax6 (Developmental Study of Hybridoma Bank, Univer- 
sity of Iowa, Iowa City, lA, USA); rabbit anti-GABA 
(Sigma); mouse anti-Pou4fl (Santa Cruz, Biotechnology, 
Dallas, TX); rabbit anti-Dbxl, a gift from Y. Nakagawa 
(University of Minnesota, Minneapolis, MN, USA) [17]; 
Alexa secondary antibodies (Invitrogen). To compare the 
spatial expression domain of molecular markers, ISH and 
immunohistochemistry were performed on serial sections 
of the same embryo or carefully matched sections of differ- 
ent embryos. 
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